AD-A141  461 
UNCLASSIFIED 


DETECTION  OF  RANDOMLY  OCCURING  SIGNALS  USING  SPECTRA  1/1 

AND  FREQUENCY  DOMAIN.. <U)  NAVAL  UNDERWATER  SYSTEMS 
CENTER  NEW  LONDON  CT  NEW  LONDON  LAB..  R  F  DWYER 
23  MAR  84  NUSC- TM- 84 1057  F/G  12/1 


NL 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS- 1963-A 


1 


TM  No.  841057 


CO 


NAVAL  UNDERWATER  SYSTEMS  CENTER 
NEW  LONDON  LABORATORY 
NEW  LONOON,  CONNECTICUT  06320 


< 

I 

D 

< 


Technical  Memorandum 


DETECTION  OF  RANDOMLY  OCCURRING  SIGNALS  USING  SPECTRA 
AND  FREQUENCY  DOMAIN  KURTOSIS  ESTIMATES 


Date:  23  March  1984 


Prepared  by:  ^ 

Roger  F.  Dwyer  ~ 


Approved  for  public  release;  distribution  unlimited 


84  05  23  017 


TM  No.  841057 


NAVAL  UNOERWATER  SYSTEMS  CENTER 
NEW  LONOON  LABORATORY 
NEW  LONOON,  CONNECTICUT  06320 


Technical  Memorandum 


DETECTION  OF  RANDOMLY  OCCURRING  SIGNALS  USING  SPECTRA 
AND  FREQUENCY  OOMAIN  KURTOSIS  ESTIMATES 


Date:  23  March  1984 


Prepared  by:  /Cocyct  j 

Roger  F.  Dwyer 


DTIC 

ELECTEJ 
MAY  2  3  1984 


Approved  for  pub’ic  release;  distribution  unlimited 


TM  No.  841057 


I 


! 


ABSTRACT 

Several  detection  statistics  are  compared  in  the  frequency  domain  based 


on  the  asymptotic  probability  of  detection  criterion.  ^fttey  include, 
secona-order ,  fourth-order,  and  two  forms  of  kurtosis  estimates.  The  results 
show  that  for  randomly  occurring  signals  or  non-Gaussian  signals,  the 
fourth-order  and  kurtosis  estimates  can  have  higher  asymptotic  probability  of 
oetection  levels  compared  with  second-order  estimates.  But,  only  for  the 
kurtosis  estimates  do  the  results  seem  significant.  Moreover,  it  a 
second-oroer  estimate  of  the  noise  is  available  to  normalize  a  fourth-order 
estimate  of  signal  and  noise,  the  resultant  moaifieo  kurtosis  estimate  has 
higher  asymptotic  probability  of  detection  levels  even  for  Gaussian  signals. 
This  result  only  holds  when  there  is  a  significant  positive  covariance  between 
the  numerator  and  the  normalizing  noise  sample  in  the  denominator.  On  the 
otner  hand,  if  an  independent  noise  sample  is  used  to  normalize  a  second-oroer 
or  fourth-order  estimate  the  overall  performance  based  on  the  asymptotic 
probability  of  detection  will  be  degraded  compared  with  the  unnormalized 
second-order  or  fourth-order  estimates,  respectively.  This  result  coulo 
impact  current  sonar  processing  methods.  v -  _ 
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I.  INTRODUCTION 

In  this  paper  the  performance  of  several  detection  statistics  are 
compared  based  on  the  asymptotic  probability  of  detection  criterion.  The 
likelihood  ratio  (LR)  is  considered  a  general  method  for  deriving  optimum 
detection  statistics  or  receiver  structures.  However,  the  particular  form  of 
the  likelihood  ratio  may  not  be  practical  to  implement,  or,  on  the  other  hand, 
the  circumstances  under  which  the  LR  is  designed  may  change.  In  these  cases 
it  may  be  of  interest  to  consider  detection  statistics  other  than  second-order 
which  may  have  better  performance  characteristics  for  non-Gaussian  signals. 

The  emphasis  is  on  detection  here,  but,  kurtosis  estimates  can  also 
distinguish  between  Gaussian  and  non-Gaussian  signals.  In  applications  this 
is  an  important  property. 

To  motivate  our  analyses  to  follow  the  form  of  the  LR  will  be  derived 
under  conditions  which  can  be  supported  by  physical  evidence.  But,  as  has 
been  suggested  already,  the  actual  LR  statistic  for  any  particular  condition 
will  not  be  analyzed  in  the  paper.  Rather,  we  will  be  interested  in  a 
statistic  which  may  be  applied  to  a  wide  range  of  situations  where 
non-Gaussian  signal  may  arise. 

Let  there  be  available  a  set  of  N  independent  and  identically  distributed 
(i.i.d.)  samples,  as  depicted  in  figure  1.  The  joint  probability  density 
function  (pdf)  for  these  samples  may  be  written  as 
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Equation  (1)  will  represent  the  noise  only  hypothesis  and  will  be  denoted  by 


Now,  if  in  the  set  of  N  independent  samples  some  belong  to  another 
probability  distribution,  say  F2U},  representing  the  signal  and  noise  and 
identified  by  H^,  as  depicted  in  figure  (2),  then  the  livelihood  ratio  for 
the  N  samples  would  be  given  by 


LR  =  g(xj,  X2,  •  •  •»  x^)/f(x^,  X£*  •  •  »*  x^)  (2) 

wnere  g(x^,  x2,  .  .  xn)  represents  the  joint  pdf  for  signal  and 
noise.  Before  giving  the  form  of  equation  (2),  consider  the  case  when  only 
one  sample  belongs  to  F2(x)  out  of  a  possible  set  N. 


In  tnis  case,  from  the  law  of  total  probability  the  joint  pdf  of  signal 
and  noise  is  given  by  for  independent  samples 


9(x^,  x2’  •  *  =  f^(x^)  f  1  ( x ^ ) 


+  P2  f^x^fgCxg)  »  f^x.) 

i  =3 


K-l 


+  PV  *  Mxi)  Mxl<)  *  Mxi) 
K  i-1  1  1  1  *  j«k+l  1  J 


+  PN  W  A  W  ’ 


(3) 
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where,  Pk>  k  =  l,  2,  .  .  N,  is  the  probability  that  the  k-th  sample 
belongs  to  Fj^).  By  forming  the  LR  between  equation  (3)  and  equation 


(1)  we  obtain 


lr  pk[f2(xk)/fi(xk)]  • 


Thererore,  the  LR  of  equation  (4j  is  formed  by  averaging  the  inaivioual 


sample  likelihooa  ratios  over  all  possible  sample  positions  of  the  signal.  If 


it  is  known  when  the  signal,  belonging  to  F2(x;,  occurs,  say  the  k-th,  then 

Pk  =  1  ana,  therefore,  LR  =  f2lxk)/fi(X|<).  Thresholding  is 

sometimes  usea  to  establish  when  the  sample  belonging  to  F2(x),  i.e.,  the 

signal,  occurs.  But  this  method  will  not  be  considerea  here  since  it  appears 


to  presuppose  the  solution  before  it  can  be  derived  from  fundamental 


principles. 


The  likelihooa  ratio  of  the  more  general  case  of  M  samples  belonging  to 


F2(x)  can  be  written  as  follows: 


N-(M-l)  N-^M-2)  N-l 

E  -E 


9  *  •••>  Xiy  )  9  (5) 


s'  ’  "K2 


Kj.I  KjJCj.l  Vl‘  Sl-2*1  Wl*1 


where,  f^l  x^  ,  x^  ,  . . . ,  x^  )  =  ( x^  j  1 2  (  x^  )  ...  ^2  ^  ^  1  ^  XK  ^  ...  ^2  ^ 


1  2 


and  it  was  assumed  that  the  probability  of  any  particular  conf iguration  was 


equally  likely,  i.e.,  v  u  =  Mi(N  -  M)I/Ni.  Therefore,  the 

l'l  *•  M 


of  equation  (b)  is  formed  for  the  N  sampling  positions  by  averaging  over  all 


possible  combinations  of  the  M  samples  containing  the  signal 
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If  M  *  1  equation  (5)  reduces  to  equation  (4)  with  =  i/n.  However, 
if  M  »  N  then  equation  (5)  reduces  to 

N 

LR  -  Cf2(x1)/f1(x1)]  ,  (6) 

which  is  the  form  of  the  likelihood  ratio  treated  by  Marcum  [1,  p.  209]. 

These  results  can  also  be  extended  to  include  dependent  samples.  For 
example,  for  Markov  dependence,  equations  (4)  and  (5)  will  be  functions  of 
conditional  probability  densities  and  the  LR  structure  will  be,  therefore, 
more  complicated. 

From  these  cases  it  appears  that  the  form  of  the  LR  depends  upon  what  is 
known  about  the  occurrence  of  the  signal  over  the  interval  in  addition  to  the 
usual  requirements  of  knowledge  of  the  pdf's.  If  we  design  our  receiver  based 
on  equation  (6)  and  the  signal  occurs  only  for  a  percentage  of  the  N  samples 
then  we  would  expect  the  performance  of  the  receiver  to  be  degraded  over  what 
could  be  achieved  under  the  proper  likelihood  ratio  formulation.  Therefore, 
the  LR  based  on  equation  (6)  may  not  lead  to  the  optimum  detection  statistic  • 
for  randomly  occurring  signals.  This  has  been  shown  to  be  the  case  by 
Ferguson  [2]  where  skew  and  kurtosis  estimates  were  optimum  statistics  under 
suitably  chosen  conditions.  Since  in  practice  the  previously  derived  LR 
structures  may  be  too  costly  to  implement  this  paper  will  compare  various 
forms  of  the  kurtosis  estimate  with  a  second-order  estimate.  These  statistics 
will  be  simpler  to  implement  but  also  may  lead  to  improved  performance. 
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II.  ASYMPTOTIC  PROBABILITY  OF  DETECTION 


The  introduction  considered  signals  which  could  occur  randomly  over  a 
sample  set  N  in  a  one-dimensional  space,  or  only  over  discrete  temporal 
locations.  Often  in  applications  it  is  desired  to  operate  on  signals  in  the 
frequency  domain.  Tnis  domain  would  allow  signals  to  occur  randomly  in 
two-dimensions,  vis.,  frequency  and  temporal  locations.  An  application  to  the 
ice-induced  signals  was  discussed  by  Dwyer  in  reference  [3].  This  idea  can 
also  De  extended  to  three-dimensions  by  taking  into  account  the  spatial 
locations.  Here  we  formulate  the  problem  in  the  frequency  domain  but  the 
basic  properties  would  be  applicable  in  any  dimension.  The  higher  dimensions, 
nowever,  would  offer  more  opportunities  for  the  signal  to  occur  randomly. 


Let  x ( i , q )  =  x{[i+(q-l)M]n}  ,  i*0,l, . . . ,M-1,  q  «  1,2, ...,n. 


represent  real  discrete  data  with  h  representing  the  interval  between 
consecutive  samples.  The  discrete  Fourier  transform  (DFT)  is  defined  as 
follows: 


x(q,F  )  =^h7R  £  W ,x(l,q)  exp(-jF  i), 

M  i-0  v 

where,  j  =*-^-1,  Fp  „  2»fpn  is  the  p-th  radian  frequency  component, 
p  *  0,l,...,M-i  and  fp  «  p/Mh  Hz.  For  simplicity,  let  the  window  weights 
equal  one,  i.e.,  *  l,  i«0,l,...,M-l,  and  h  =»  1. 


<  •  if  s  l-fr. : 
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The  time  domain  data  will  be  represented  as  an  additive  sum  of  signal  and 
noise  of  the  form, 

x(i,q)  =  N(i,q)  +  m(i,q)  S(i,q)  (8) 

where  the  noise  N(i,q)  and  signal  S(i,q)  are  zero-mean  Gaussian  process.  The 
function  m(i,q)  modulates  the  signal  in  such  a  way  that  the  frequency  domain 
representation  will  be  a  non-Gauss ian  process  [3]. 

Nuttall  [4]  considered  a  signal  that  was  modulated  by  a  random  constant 
and,  tnerefore,  was  not  a  function  of  time.  This  kind  of  modulation  as  was 
found  by  Nuttall  would  not  lead  to  a  non-Gaussian  process  based  on  the  model 
and,  therefore,  the  LR  of  equation  (6)  would  be  appropiate.  In  contrast,  the 
modulation  function  of  this  paper  is  a  random  function  of  time  over  the 
detection  interval  and  leads  to  non-Gaussian  process.  This  model  is  supported 
by  real  data  measurements  as  discussed  in  reference  [3]. 

The  power  spectrum  density  (PSD)  which  represents  a  second-order  estimate 
is  defined  as  [5,6] 

P(FJ  -  (1/n)  £  x(q,F  )x*(q,F  ) ,  (9) 

w  q«l  K  * 

where  the  asterisk  represents  complex  conjugate.  The  asymptotic  variance  of 
power  spectra  and  complex  cross  spectra  were  discussed  in  reference  [7].  The 
parameter  n  is  sometimes  called  the  degrees  of  freedom  of  the  estimate. 
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Let  the  frequency  domain  Kurtosis  (FDK)  estimate  [3]  be  defined  for  real 
and  imaginary  parts  of  equation  (7)  separately.  We  can  therefore  discuss  the 
tneoretical  properties  of  the  FDK  for  its  real  part  only.  The  imaginery  part 
has,  identical  properties.  In  appl  ications  both  real  and  imaginary  parts  are 
estimated,  since  both  contain  information.  In  addition,  the  real  and 
imaginary  parts  of  the  FDK  estimate  can  be  combined  to  form  one  statistic, 
say,  for  example,  the  magnitude.  But  these  considerations  are  not  expected  to 
change  tne  conclusions  to  follow  in  any  significant  way.  Although,  a 
performance  improvement  based  on  the  magnitude  of  the  FDK  estimates  over  the 
real  or  imaginary  part  separately  would  be  expected  compared  with  the  PSD.  An 
example  for  the  magnitude  of  the  FDK  estimates  will  be  given  later. 

Tne  FDK  estimate  for  the  real  part  of  equation  (7)  is  defined  assuming 
tnat  x(q,Fp)  is  a  zero-mean  process  as 

"  A  >2 

K(FJ  -  (1/n)  £  [x(q,  F  ) J4/  (1/n)  £  [x(q,F  )]2  .  (10) 

v  q=l  v  q=l  K  ) 

If  x(q,Fp)  is  not  a  zero-mean  process  then  the  mean  would  be  estimated 
from  the  data  and  subtracted  from  x(q,F  ).  So,  therefore,  the  mean  can  be 
accounted  for  in  principle  and  need  not  concern  us  further. 

The  asymptotic  probability  of  detection  (APD)  introduced  in  reference  [8] 
will  oe  deri^d  for  three  specific  cases  of  equation  (10)  and  compared  with 

the  APD  results  of  equation  (9).  In  case  1,  the  numerator  and  denominator  of 
equation  (10)  will  contain  signal  and  noise  as  shown  in  the  data  model  of 
equation  (8).  This  was  the  form  of  the  frequency  domain  kurtosis  estimate 
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used  in  tne  real  data  analysis  of  reference  [9J.  For  case  2,  the  denominator 
will  contain  a  noise  only  estimate.  The  results  of  this  case  are  of 
tneoretical  interest.  However,  some  comments  concerning  implementing  this 
case  will  be  discussed  later.  In  the  last  case  the  denominator  of  equation 
(10)  will  be  set  equal  to  one.  Thus,  case  3  represents  a  fourth-order  moment 
estimate. 


Since  only  tne  real  part  will  be  treated  theoretically,  we  shall  write 
equation  (7)  as  follows: 


xlq.Fp) 


(1/M)i/2 


M-l 

Y  N(i,q)  cos  (F  i) 
i=0  p 


M-l 

m(q)  Y  S(i,q)  cos  (F  i)  , 
i=Q  p 


n 


where  m(i,q)  was  assumed  to  change  slowly  with  i  and  therefore,  could  be  w  , 
approximated  as  m(i,q)  =  m(q). 


To  simplify  notation  equation  (11)  will  be  expressed  for  a  particular 
frequency  as. 


x(q)  -  N(q)  +  m(q)S(q) .  (12} 

where  N(q),  m(q),  and  S(q)  are  mutually  independent.  In  addition,  x(q)  will 
be  assumed  a  zero-mean  process  and  statistically  independent,  i.e., 

E[x(q1)x(q2)]  -  «(qi-q2)  £[x(q1)2].  Where  ^(qpqj)  is  the 
Kronecker  delta  function.  The  independent  assumption  is  needed  in  order  to 
evaluate  the  variance  of  equations  (9)  and  (10).  However,  the  time  series 
x ( 1 ,q )  with  respect  to  i  may  be  considered  dependent.  This  point  will  be 
discussed  more  fully  later. 
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Tne  modulating  term  m(q)  will  be  modeled  as  a  oernoulli  time  series  [3] 
defined  as  follows: 

1  ;  Pr{>n(q)  =  1}  =  L 

m(q)  = 

0  ;  Pr{m(q)  =  0  }  =  1  -  L. 

Therefore,  all  moments  of  m(q)  are  given  by  E[m(q)r]  =  L, 
r  -  1,  2,  ...  .  Thus,  m(q)  modulates  the  signal,  in  the  frequency  domain  at  a 
specific  frequency,  oy  turning  it  on  or  off  over  the  detection  interval  .  The 
probability  (L)  of  the  signal  Deing  on  will  be  a  parameter  of  the  performance 
resul ts. 

The  asymptotic  (n  *■«)  behavior  of  equation  (9)  and  (10)  can  be  obtained 
from  Cramer's  convergence  proofs  [10]. 

For  a  function  P(x^,x2z),  z  =  1,  2,  3,  of  two  central  moments  x^, 
x2z’  ^rotn  a  one-dimensiona1  sample  corresponding  to  equation  (9)  and  (10) 
Cramer  has  shown  tnat  as  n  *«,  F  »  n{E(F) ,var(F) f  ,  where  n{  ,  } represents  a 
normal  process  with  mean  and  variance  given  by,  respectively, 

E[F]  =  F[E(x1),£(x2z)3  +  0(l/n), 

VAR[F]  .  VAR(Xl)F*  +  VAR(x2  ) 

*  2F1F2C0V(x1  x21)  +  0(1 In312). 

The  parameters  F^  and  F2  are  partial  derivatives  of  F,  i.e.. 
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In  tne  following  analysis  we  shall  replace  F  by  P  for  the  PSD  and  replace 
F  by  K  for  the  FDK  when  evaluating  their  respective  means  and  variances. 

Since  tne  asymptotic  process  is  Gaussian,  the  asymptotic  probability  of 
detection  (APD)  can  be  written  as  follows: 

APO  ~  1  -*j*-1(l-a)  a0(F)/Ol(F)  +  £EQ(F)  -  EjfFJj/a^F) }  ,  (13) 

where,  E..{  ),  a^(  )  are,  respectively,  the  mean  and  standard  deviation  of 
the  noise  only  process  when  i  =0,  or  the  signal  and  noise  process  when 
i  =  1.  [is  tne  standard  error  function  and  a  represents  the  desired 
fal  se  al  arm  probabil  ity. 

III.  THEORETICAL  DETECTION  PERFORMANCE  RESULTS 

A.  Power  Spectrum  Density 

Tne  power  spectrum  density  (PSD)  estimate  for  the  real  part  of  equation 
(9)  is  given  by 


P(F  )  =  F(xx,l)  *  (1/n)  Z  *(q)  •  (14) 

M  q=l 

The  required  components  of  tne  APD  are 
ELP(Fp)3  -  N(1  +  L  SNR),  VAR  [P(Fp)]  =  (1/n)  2N2  a^ 
where,  N  *  E[N(q)2],  S  =  E[S(q)2],  SNR  =  S/N,  and 
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Recall  that  the  variance  of  the  periodogram,  as  well  as  equation  (14), 
does  not  decrease  with  increasing  transform  size  M  [6].  We  shall  prove  later 
that  the  variance  of  the  FDK  is  also  independent  of  M  and  like  the  PSD  both 
variances  decrease  proportional  to  (1/n)  if  n  sufficiently  large. 

By  substituting  the  above  parameters  for  the  PSD  into  equation  (13)  we 
obtain 


APD  =  1  -  ^<*>Qt>-1(l  -  a)  -  L  SNRv'"n/vr2]/v/T1| .  (15) 

Equation  (15)  will  be  evaluated  for  specific  values  of  the  parameters  and 
compared  with  the  results  for  the  FDK  of  case  1  in  the  next  section. 

Another  point  should  be  mentioned  concerning  the  PSD  before  proceeding. 
Later  we  shall  discuss  the  performance  of  a  fourth  order  estimate  normalized 
by  a  noise  only  second-order  estimate,  which  is  designated  as  case  2.  This 
procedure  leads  to  higher  APD  values  compared  with  all  other  cases  including 
the  PSD  even  for  Gaussian  processes.  It  is  believed  the  reason  for  this 
depends  on  the  noise  only  second-order  estimate  normalizing  function.  But, 
the  PSD  can  also  be  normalized  by  a  noise  only  second-order  estimate.  When 
this  is  done  it  appears  that  the  PSD  would  have  a  larger  APD  for  Gaussian 
proceses.  However,  the  false  alarm  probability  cannot  be  controlled  at  a 
satisfactory  small  level  but  is  fixed  at  .5.  A  proof  of  this  is  given  in 
Appendix  A. 
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B.  Frequency  Domain  Kurtosis,  Case  1 

From  the  results  of  the  previous  section  the  FDK  can  be  expressed  as 
foil ows: 


W  =  *l'x2z  * 


(1 


where. 


n 

x,  =•  (1/n)  ^  x(q)4  , 
1  q=l 


2 

X2z  *  C(l/n) 

q*l 

Therefore,  for  z  ■  1,  case  1,  the  asymptotic  expected  value  of  equation 
(16)  reduces  to 


lim  ECkj^F  )]  -  E [x^x^J  =  E[x(q)4]/E[x(q)2]  . 
n*  oo  ^ 

Since  all  our  results  are  asymptotic  (n  **>)  we  shall  not  explicitly 
indicate  the  limiting  process  in  the  following  expressions. 

Using  tne  data  model  of  equation  (12)  the  expected  value  of  k^(Fp) 
can  De  put  in  the  parmeterized  form 

E[k,  (F  )]  =*  3  a?/a2. 
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where,  a2  -  1  +  2L  SNR  ♦  L  SNR2,  and  a3  .  1  +  L  SNR. 

The  derivaton  of  the  variance  of  equation  (16)  requires  evaluating 
several  components.  Some  of  them  are  tedious  calculations.  Therefore,  the 
main  points  of  the  derivation  will  only  be  outlined  here.  The  details  of  the 
derivation  will  be  found  in  Appendix  B. 


Based  on  Cramer's  [10]  convergence  proof  the  form  of  the  variance  of 
^(Fp)  is  given  by 


VARClc1(Fp)]  -  VAR  (xL)  F2  ♦  VAR(x21)  F2 


+  2  Fx  F2  cov(x1  x21)  , 


where, 


F1  3  (3  kl(Fp>/axl> 


X1  "  E<xl> 


E[x(q)2]"2  , 


Fo  3  )/3  ) 

2  1  ’  21  |x2l  -  E(,a) 

The  other  components  are  given  by 


-E[x(q)4]  E[x(q)2]'  . 


VAR(xx)  =  (1/n)  VAR[x(q)4] 

VAR(x21)  -  Ej[(l/ri)  £  x(q)2]4!  -  j  Ejc(l/n)  £  x(q)2]2} 
'  q=l  \  1  q=l  j 

cov[x1  x21]  -  E(xx  x21)  -  E[xl3  E[x21J 

n  n  n 


E  j(l/n3)  £  2  2  x(qx)4  x(q2)2  x(q3)2J 
«  q^«l  q2«l  q3»l  ) 
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n 

-ELx(q)4]  Ej[(l/n)]T  x(q)2]2  J 
q-1 


By  expanding  the  components  further  we  obtain. 


VARO^Fp)] 


+  4 
-  4 


-  (1/n) 


VAR[x(q)4]/E[x(q)2]4 


E[x(q)4]  E[x(q)2]2  -  E[x(q)2]4  !  E[x(q)4]2/E[x(q)2]8 

ew<,,6]  Ew<,)2]-Ecx(”4]E[^a,2]2l,E[’<i’,2]i 


Substituting  the  data  model  of  equation  (12)  into  the  above  expressions 
the  desired  parameterized  form  for  the  variance  of  case  1  reduces  to  the 
f ol  1 owi ng . 


VARCk^Fp)]  =  (1/n)  Qx  , 


(17) 


wnere. 


Qj  *  96  +  72  02  -  144  03 

^1  *  ®4^a2 

a4  >  1  ♦  4L  SNR  +  [612L  -  36L2)/96]  SNR2  +  [(420L  -  36L2)/96]SNR3 
+  C(10SL  -  9L2)/96]  SNR4 

02  -  (a|/aj)  (1/8)  [18  -  10  a3  a5/a2] 
a5  -  1  ♦  2L  SNR  ♦  [4L2  ♦  6L)/10]  SNR2 
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03  -  (a|)a|)(l/12)  [15  a6(l  ♦  L  SNK)/a§  -  3  a3/a2J 

ag  »  1  ♦  3L  SNR  +  3L  SNR2  +  L  SNR3. 

Under  the  noise  only  case  (SNR  =  0)  or  when  L  equals  0  or  1  the  variance 
of  case  1  asymptotically  approaches, 

VARCk^FpJ/Ho]  -  24/n  (18) 

Tnis  is  the  result  ootained  for  a  Gaussian  process  by  Pearson  [11]. 

Now  we  are  in  a  position  to  express  the  APD  for  the  FDK  of  case  1.  The 
result  is  given  by 

APO  -  1  -  *{^4  ♦“i  (1  -  a)  +  3(1  -  a2/a3)yF]/Q{/2f.  (19) 

For  L  equal  to  0  or  1  the  APO  of  equation  (19)  equals  the  false  alarm 
probability  so  detection  is  not  possible  under  these  conditions.  Detection  is 
also  not  possible  for  the  PSD  when  L  equals  zero.  But  unlike  the  FDK  of  case 
1  tne  probability  of  detection  for  the  PSD  is  maximum  when  L  equals  one.  The 
FDK  of  case  2  will  also  have  this  property.  This  means  that  x^Fp)  as  a 
detection  statistic  is  only  sensitive  to  non-Gaussian  (0  <  L  <  1)  signals. 

An  example  will  demonstrate  these  results  more  clearly.  Figure  3 
represents  a  comparison  of  the  PSD  and  FDK  of  case  1  as  a  function  of  the 
probability  of  occurrence  (L)  of  the  signal  for  several  SNR  values.  The 
figure  snows  that  the  FDK  can  nave  a  higher  probability  of  detection 
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depending  on  SNR  for  small  values  of  L  only.  As  L  increases  the  probability 
of  detection  for  the  PSD  increases  with  L  and  will  surpass  the  probability  of 
detection  for  the  FDK  of  case  1.  As  L  approaches  one  the  probability  of 
detection  for  the  FDK  of  case  1  decreases  and  approaches  the  false  alarm 
proDability  (a). 

A  real  data  example  will  put  these  results  in  perspective.  The  following 
example  has  been  explained  in  reference  [3].  It  is  from  Arctic  under-ice 
environmental  data.  Figure  4  compares  the  power  spectrum  density  estimate  and 
the  real  part  of  the  frequency  domain  kurtosis  estimate  for  a  segment  of 
Arctic  data.  Both  estimates  are  averaged  over  many  consecutive  FFT 
estimates.  We  shall  concentrate  on  the  frequency  with  the  highest  kurtosis 
estimate  and  a  corresponding  small  PSD  estimate  on  the  figure.  From  the 
theoretical  results  we  know  that  the  signal  causing  the  high  kurtosis  level 
only  occurs  a  small  percentage  of  time  and  when  it  occurs  it  probably  has  a 
high  SNR  level.  This  conclusion  could  not  be  reached  from  the  PSD  estimate 
alone.  Therefore,  the  FDK  estimate  of  case  1  contributes  additional 
information.  So,  the  usefulness  of  the  FDK  estimate  of  case  1  cannot  be 
determined  from  its  probability  of  detection  performance  level  alone.  But 
must  be  evaluated  in  terms  of  contributing  additional  information  which  the 
PSD  estimate  is  insensitive  to. 

For  the  next  two  cases  of  the  FDK  we  will  explore  ways  of  improving  its 
probability  of  detection.  But,  on  the  other  nand,  this  procedure  may 
sacrifice  the  information  content  of  the  FDK  estimate. 
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C.  Frequency  Domain  Kurtosis,  Case  2 

For  this  case  the  denominator  of  equation  (10)  contains  a  noise  only 
second-order  estimate.  The  FDK  estimate  of  case  2  follows  from  equation  (16) 
as 


k2<Fp)  a  xl/x22.  (20) 

The  derivations  of  the  following  expressions  are  similar  to  the 
derivations  for  the  FDK  of  case  1.  Therefore,  only  the  main  points  of  the 
derivation  will  be  explicitly  stated  in  the  following. 

For  sufficiently  large  n,  the  expected  value  of  equation  (20)  is  given  by 

E|>2(Fp)]  -  E[x(q)4/H1]/  £[x(q)2/H0]  2  (21) 

Where  we  now  signify  which  hypothesis  the  expected  value  is  conditioned 
on.  8y  substituting  tne  data  model  of  equation  (2)  into  equation  (21)  we 
obtain  the  parameterized  result. 


El>2(Fp)]  -  3  a2  (22) 

The  variance  is  given  by 

VAR|>2(Fp)]  -  (l/n)|  VAR[x(q)4/H1]/E[x(q)2/H0]2 

+  4  E[x(q)4/H1J2|E[x(q)4/H0]  E[x(q)2/HQ]2  -  E[x(q)2/H0]4f/E[x(q)2/H0] 
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~2{E[x(q)4/H1]  /  E[x(q)2/H0]6}j 

n  n 

E[(l/n)  [x(q)4/H1]  [(1/n)  £  Cx(q)2/HQ]]2 

q-1  q=l 

n  n 

-E[(l/n)  2  x(q)4/Hi]  E[(l/n)  ^  x(q)2/HQ]2 
q»l  q=l 

The  last  term  of  equation  (23)  represents  tne  covariance  of  x^  and 
X£2»  i.e,  cov[x]X22]-  If  tne  noise  only  estimate  X22  was 
statistically  independent  of  x^  tnen  the  covariance  would  equal  zero.  This 
may  happen  if  tne  noise  estimate  is  from  an  adjacent  frequency  location  or 
from  another  interval  of  time.  This  would  change  the  APO  performance  results 
consideraoly. 

Employing  tne  data  model  we  obtain  the  parameterized  result 


VAR[k2(Fp)]  =  (l/n)  Q2 


(24) 


where  Q2  -  96  a4  +  72  a2  -  144  a2  a3. 

If  tne  covariance  is  zero  tne  tnird  term  of  Qj  would  be  zero  giving  a 
higher  variance  for  Now  the  role  of  the  normalizing 

component,  Xj2,  is  clear.  If  X22  is  correlated  with  the  noise  component 
of  x1  tnen  tne  effect  is  to  reduce  the  variance  of  ^2 ( Fp ) • 

Tne  APD  of  case  2  also  shows  this  relationship. 

APD  -  1  (1  -  a)  -  3  a 


(25) 
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where  a?  -  2L  SNR  ♦  L  SNR2. 

Therefore,  if  Q2  is  minimized  the  probability  of  detection  can  be 
improved.  Tne  APO  is  maximum  for  the  FDK  of  case  2.  In  practice  a  partially 
correlated  noise  sample  might  be  available  to  improve  the  detection 
performance  of  the  FDK  estimate.  An  example  will  be  given  later  to 
demonstrate  the  idea.  But  before  given  examples  of  the  results  of  case  3  will 
be  given. 

D.  Frequency  Domain  Kurtosis,  Case  3 

Recall,  for  tnis  case  the  denominator  of  equation  (10)  is  set  equal  to 
one.  Tne  derivations  of  the  mean  and  variance  of  case  3  are  similar  to  the 
previous  cases  so  we  will  only  present  the  results. 

The  expected  value  and  variance  for  the  FDK  of  case  3  are  given  by 

E[K3(Fp)]  =  3  N2  a2  , 

VAR[K3(Fp)]  =  ( 1 / n ) ^96  N4  a4) 

Tne  APD  follows  as 

APD  =  1  -*j[<r1(l  -  a)  -  3  a7 /r/v56]//a“  J  .  (26) 

At  tnis  point  it  is  appropriate  to  consider  the  performance  of  a  V-th 
order  power  law  using  tne  data  model  of  equation  (12).  Recall  that  both  N(q) 
and  S(q)  are  assumed  to  be  statistical  independent  and  zero-mean  Gaussian 
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processes.  Here  we  shall  utilize  one  degree  of  freedom,  so  the  results  are 
not  asymptotic. 


The  V-th  order  power  law  is  defined  as  follows 


y  =  xv 


where  for  v  even,  v  »  2k,  k  =  1,2,...,,  and  for  v  odd,  v  =  2k— 1 ,  k-1,2,...  . 


When  v  equal s  two  the  results  will  apply  to  equation  (14)  and  when  v 
equals  four  they  will  apply  to  equation  (16),  where  X£3  „  i. 


Case  (I)  V  EVEN 


By  rewriting  equation  (27)  in  the  following  way 


y  exp(jn2«)  =  xv,  n»0,+l,+2, . . . , 


the  solutions  for  x  can  be  obtained.  For  y  >  0 


x  »  yl/v  exp(jn2ir/v) . 


By  examining  the  possible  solutions  we  find  that  there  are  only  two  real 
solutions  for  x,  viz.. 


*  y^v  and  xg  -  -y* 


-vriinritfri  ir 
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There  are  no  real  solutions  for  x  when  y  <  0.  Therefore,  the  probability 
density  function  of  y  is  given  by 

Vy)  =  Lvy{v_1)/vr1  Cfx(y1/V)  ♦  fx(-y1/v)J,  y  >  0 
fy(y7  =  0  ,  y  <  0. 

From  the  characteristic  function  of  equation  (12)  the  input  probability 
density  function  can  be  shown  to  be  given  by  the  following  mixture  density 


fx(x)  =  ( 1— L )  fx(x)  ♦  Lf2(x), 


(28) 


where 


,  f^x)  is  n|o,  on|,  and  f2(x)  is  n|o,  oj  ,  a  =  +  “s^2. 


Therefore,  since  a  Gaussian  process  is  symmetric  about  its  mean  value  we 
obtain 

My)  -  2Lvy^-1)/vr1  ju-l)  fi(y1/v)  +  Lf2  (y1/v)j  ,  y  >  0 


and 


Fy(y)  »  2|(l-L)L4>(y1/v/on)  -  1/2]  +  L[<J>(y1/v/o)  -  l/2]j,  y  >  0 


where,  <t>(  )  is  the  Gaussian  error  function. 
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Since  we  are  interested  in  the  performance  properties  with  one  degree  of 
freedom,  we  shall  let  L=l,  i.e.,  the  signal  occurs  with  probability  one  on 
this  sample.  The  objective  is  to  find  the  probability  of  detecting  this 
signal  at  the  output  of  a  v-th  order  power  law  for  a  fixed  false  alarm 
probabil  ity. 

Tne  false  alarm  probability  is  defined  as 

3  -  1  “  FY<ya/H0),  y  >  0 

where  y^  is  the  threshold  (value  of  y)  where  the  desired  false  alarm 
probability  is  maintained.  For  tnis  case  we  find 

ya/V  =  ®n  -  a/23- 

The  probability  of  detection  (Pd)  is  therefore 

Pd  =  1  -  Fy(ya/H1) 

=  2jl  -*[+“1(l-a/Z)/(l+SNR )1/2](,  y  >  0  (29) 


wnere  SNK  ■  o2/02. 

Once  a  is  fixed  for  each  v  the  probability  of  detection  is  independent  of 
v  and  only  depends  on  SNR.  Therefore,  the  performance  of  all  even  order  power 
laws  will  be  the  same  assuming  Gaussian  processes  based  on  the  probability  of 
detection  criterion. 
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Case  (II)  V  ODD 

The  odd  order  power  law  is  given  by 

y  =  xv,  v  =  2k-l,  k=l,2,...  . 

For  y  >  0  and  x  >  0  the  solutions  of  x  are  represented  by 
x  =  yl/v  exp(  jn2if/v) . 

But  there  is  only  one  real  solution,  viz.,  x  =  y^v. 


For  x  <  0  and  y  <  0,  the  odd  order  power  law  is  given  by 

-|y|  *  -|x|v»  or  y^  =  xv:  under  the  conditions 

0<y^<oo,  0<X|<oo 

-<» <  y  <  0,  -»<  x  <  0 

i 

So,  the  probability  density  function  at  the  output  is  given  by 
fy(y)  -  [V  y(,-u/v]'1  fx(yl/v),  y  >  0 

'„<*!»  -  C.  y1(,'1)/v]"1  y,  >  0 

Integrating  we  obtain 


Fy(y) 


fy(Z)dZ  = 


/  yz)d‘ 


fyU)dZ 
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1/2  + 


/  (2ir )  1 

•7  n 


expC-g^/a)  dy,  or 


Fy(y)  =*[y1/v/a],  y  >  o  . 

If  we  fix  a  threshold  to  the  right  (y  >  o)  we  find  for  «, 


«  “  1  -  Fy(ya/H0)  and 


$  ^(l-«) 


Under  this  threshold  the  probability  of  detection  is 


Pd  =  1  -  Fy(ya/n1)  =  1  -*[*-1(l-a)A/l+SNR] 

which  is  also  independent  of  v.  As  SNR  approaches  infinity,  for  ya  finite, 

we  obtain,  1  im  Pd  >  1/2. 

SNR-oo 

With  one  threshold  this  is  the  best  we  can  do.  But  for  two-sided  distributions 
even-order  power  laws  are  usually  employed.  However,  if  we  set  two  threshol  ds 
located  at  y  and  -y  for  y  >0,  then 

a  a  <x 

a  =  2[1  -  Fy(ya/Hg)],  which  gives  the  threshold  value  at, 
y<*  “  °n  ”  °72). 


Then 
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which  is  identical  to  the  results  for  the  even  order  power  laws. 


We  conclude  that  for  zero  mean  Gaussian  processes  there  is  no  essential 
difference  in  performance  for  the  v-th  order  power  laws  based  on  the 
probability  of  detection  criterion,  if  two  thresholds  are  set  for  the  odd 
orders.  These  results  also  hold  true  if  L  ^  1. 

An  interesting  question  presents  itself.  Can  performance  improvements  be 
obtained  by  summing  over  a  large  number  of  samples  (n  large)  compared  with  the 
single  degree  of  freedom  case?  For  randomly  occurring  signals  the  answer  is 
not  obvious,  as  would  be  so  for  purely  Gaussian  signals.  By  comparing 
equations  (29)  and  (19)  under  the  same  SNR  and  false  alarm  probability  we  find 
that  the  FDK  of  case  1  can  have  higher  asymptotic  probability  of  detection 
levels  compared  with  the  probability  of  detection  for  the  single  degree  of 
freedom  case  of  equation  (19),  depending  on  the  probability  of  occurrence  (L) 
of  the  signal.  This  conclusion  also  holds  for  the  FDK  of  cases  2  and  3. 

The  probability  of  detection  cannot  be  obtained  for  the  FDK  of  case  2 
with  one  degree  of  freedom  since  the  joint  distribution  of  the  numerator  and 
denominator  is  unknown.  This  is  also  true  for  the  FDK  of  case  1  if  the 
degrees  of  freedom  are  greater  than  two  but  not  large  enough  to  assume 
Gaussianity.  On  the  other  hand,  if  cov(x1x22)  equals  zero  than  with  only 
a  limited  number  of  degrees  of  freedom  the  probability  of  detection  can  be 
obtained. 
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We  will  now  give  several  examples  comparing  the  FDK  of  cases  2  and  3  with 
the  PSD.  Figures  5  tnrough  9  represent  a  comparison  of  the  PSD  with  cases  2 

and  3  based  on  the  asymptotic  probability  of  detection  vs  L  for  fixed  SNR  or 

*3 

vs  SNR  for  fixed  L.  In  all  the  figures  the  false  alarm  probability  is  10  J 
and  the  sample  size  (n)  is  2000.  The  results  show  that  the  FDK  of  case  2  has 
a  higher  asymptotic  probability  of  detection  in  all  the  figures  even  for  a 
Gaussian  process,  i.e.,  L  -  1.  The  last  section  gave  reasons  for  this  result 
and  other  comments  are  given  in  appendix  A.  Under  some  situations  the  FDK  of 
case  3  has  slightly  nigher  asymptotic  probability  of  detection  compared  with 
tne  PSD  but  only  for  non-Gaussian  processes.  Figures  8  and  9  clearly  show  how 
the  APD  changes  for  the  FDK  of  case  3  and  the  PSD  as  L  is  varied  from  .1  to 

1.0.  For  a  non-Gaussian  process  (L  =  .1)  the  FDK  of  case  3  has  higher  APD 

levels  vs  SNR  compared  witn  the  PSD.  Whereas,  the  reverse  is  true  for  a 
Gaussian  process  (L  =  1)  as  shown  in  figure  9.  These  conclusions  do  not 
change  with  changes  in  the  sample  size  or  false  alarm  probability. 

In  practice  the  asymptotic  probability  of  detection  for  the  FDK  will  be 
closer  to  the  results  of  case  2  if  the  correlation  between  x^  and  x22  is 
hign.  On  the  other  hand,  if  the  correlation  is  zero  or  low  then  the  FDK  of 
case  3  should  be  utilized  instead  of  normalizing  the  fourth-order  estimate  by 
an  uncorrelated  noise  sample.  This  is  true  because  the  uncorrel  ated  noise 
sample  tends  to  increase  the  variance  and,  therefore,  reduce  the  APD. 

IV.  SIMULATIONS 

Several  simulations  of  the  ideas  presented  in  the  previous  sections  will 
be  given  to  check  the  theoretical  results.  We  will  utilize  equation  (12)  as 
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the  data  model  hut  here  both  real  and  imaginary  parts  will  be  included  in  the 
simulations.  Therefore,  the  real  and  imaginary  parts,  respectively,  are  as 
follows: 


XR(q)  =■  Nr(1)  +  mR(q)SR(q) 

(27) 

X^q)  *  Nj(q)  +  mjfaJSjtq) 

The  SNR  is  defined  as,  SNR  =  E[SR(q)2]/E[NR(q)2]  =  E[SI(q)2]/E[NI(q)2], 
since  botn  real  and  imaginery  parts  will  have  identical  statistics  in  the 

following  simulations.  The  power  spectrum  density  is,  therefore, 

n 

p(pp)  =  (l/n)^CXR(q)2  +  Xj  (q  )2]  ,  (28) 

q=l 

and  the  corresponding  frequency  domain  Icurtosis  estimate  will  be  defined  in 
the  following  way. 


MFp)  »  CK^(Fp)2  ♦  K*(Fp)2]  1/2  .  (29) 

Other  estimates  could  be  defined  for  the  FDK  depending  on  the  intended 
application.  But  equation  (29)  will  be  our  basis  for  comparison  with  equation 
(28)  in  tne  simulations. 

A.  Simulation  of  Case  1 

Tne  real  data  example  of  figure  (4)  has  already  demonstrated  the 
usefulness  of  the  FSK  for  case  1.  Therefore,  only  limited  examples  will  be 
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presented.  The  parameters  needed  to  define  a  particular  simulation  will  be 
denoted  in  the  form  (z,  L,  SNR,  n),  where  z  equals  1  or  2  depending  on  whether 
tne  FDK  if  case  1  or  case  2  is  being  considered.  The  FDK  of  case  3  will  not 
be  simulated  here. 

Consider  the  simulation  (1,  .005,  16.0,  200)  shown  in  figure  10.  The 
data  were  generated  using  a  1024  point  FFT.  Tne  top  graph  represents  the  PSD 
expressed  in  equation  (28).  Tne  lower  graph  is  the  corresponding  FDK  estimate 
of  equation  (29).  The  signal  only  occurs  once  for  each  of  200  frequency 
locations  starting  at  frequency  100  during  the  200  consecutive  FFT  data 
samples.  But  the  temporal  location  is  random  in  each  frequency  and  therefore 
unknown.  The  resultant  estimates  have  a  form  of  a  broadband  signal.  From  the 
raw  data  in  the  figure  we  see  that  the  FDK  estimate  identifies  the  frequency 
location  of  the  randomly  occurring  signal  but  this  information  is  not  present 
in  the  PSD  estimate.  This  result  was  predicted  in  figure  (3).  So,  the 
simulation  corroborates  the  theoretical  results  for  the  FDK  of  case  1. 

Another  possible  application  of  these  resul ts  is  the  identif ication  of 
the  track  of  a  signal  tnat  is  changing  its  frequency  with  time.  The 
instantaneous  frequency  location  could  not  be  identified  from  the  FDK  but  it 
could  possibly  be  deduced  from  the  results.  The  usefulness  of  the  FDK  in  the 
application  would  nave  to  be  carefully  compared  with  tracking  methods.  Since 
tnis  would  require  a  detailed  study  it  has  not  been  treated.  But  it  does 
appear  that  the  FDK  would  be  easier  to  implement.  So  a  study  is  probably 
warranted  in  the  future.  These  considerations  would  also  apply  for  spatially 
tracked  signals. 
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B.  Simulation  of  Case  2 

In  tne  following  simulations  we  shall  perform  a  detection  experiment. 

For  noise  only  a  threshold  will  be  founa  for  the  PSD  and  FDK  by  picking  the 
highest  level  in  each  case  of  the  512  frequency  locations.  Then  for  signal 
and  noise  present  the  number  of  frequency  locations  in  each  case  exceeding 
their  corresponding  thresholds  will  be  counted.  Starting  at  location  100 
there  will  be  100  locations  with  a  randomly  occurring  signal  present  in  this 
simulation.  Since  the  FDK  of  case  2  requires  a  noise-only  second-order 
estimate  to  normalize  the  fourth-order  estimate  of  signal  and  noise,  these 
results  may  be  of  theoretical  interests  rather  than  having  practical  value. 

On  tne  other  hand,  if  a  noise-only  second-order  estimate  is  available  which 
was  highly  correlated  witn  the  noise  in  the  fourth-order  estimate  of  signal 
and  noise  then  these  results  would  be  of  practical  value.  It  should  also  be 
pointed  out  tnat  only  a  noise-only  second-order  estimate  is  needed  ana  not 
knowledge  of  any  particular  noise  sample. 

Figures  11,  12,  and  13  represent  the  results  of  a  simulation  of  (2,  .02, 
4,200),  (2,  .04,  4,200),  and  (2,  .08,  4,200).  As  before  the  PSD  estimate  is 
in  the  top  grapn  and  its  corresponding  FDK  estimate  is  in  the  lower  graph. 

The  number  of  locations  exceeding  the  threshold  are  1,  6,  and  34  for  the  PSD 
and  47,88,  and  100  exceeded  the  threshold  for  the  FDK  in  figures  11,  12,  and 
13,  respectively. 

The  theory  also  predicts  that  the  FDK  of  case  2  will  have  a  higher 
probability  of  detection  compared  with  the  PSD  estimate  even  for  Gaussian 
processes.  This  result  was  checxed  in  the  following  simulation.  The  results 
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of  the  simulation  are  shown  in  figure  14,  unoer  the  conditions  (2,  1,  SNR, 

200)  where  SNR  equals  .1,  .16,  .25,  and  .4.  The  respective  detection  results 
are  6,  13,  45,  ana  97  for  the  PSD  estimate,  and  52,  69,  99,  and  100  for  the 
FDK  estimate.  These  results  obviously  support  the  theoretical  predictions  for 
knowleage  of  a  noise-only  second-order  estimate.  Appendix  A  points  out  that 
if  the  PSD  estimate  is  also  normalized  by  a  noise-only  secona-order  estimate 
its  probability  of  detection  would  appear  to  be  better  than  the  FDK's  but  the 
false  alarm  rate  cannot  be  maintained  at  a  desired  low  level. 

V.  SUMMARY 

Since  signals  often  appear  to  occur  randomly,  especially  in  unaerwater 
acoustic  detection  problems,  ana  the  corresponding  detector  based  on  the 
likelihooa  ratio  has  a  complex  structure,  fourtn-oraer  ana  kurtosis  estimates 
were  considered  as  alternative  processing  metnods.  These  higher-order  methods 
were  compared  witn  second-order  estimates  for  signals  which  occur  randomly  in 
time  and  frequency  and  which  could  be  described  as  a  non-Gaussian  process.  It 
was  initially  believed  that  a  secona-order  estimate  was  optimum  for  Guassian 
processes.  However,  it  was  shown  theoretically,  based  on  the  asymptotic 
probability  of  detection  that  if  a  noise-only  second-order  estimate  was 
available  to  normalize  a  fourth-order  estimate  of  signal  and  noise  then  the 
resultant  modified  kurtosis  estimate  had  higher  probability  of  detection 
levels  even  for  Gaussian  processes.  On  the  other  hand,  if  an  independent 
noise  sample  is  used  to  normalize  a  secono-oraer  or  fourth-oroer  estimate  the 
overall  performance  based  on  the  asymptotic  probability  of  detection  will  be 
degraded  compared  with  unnormalized  secona-oraer  or  fourth-oraer  estimates, 
respectively.  Only  it  there  is  positive  covariance  between  the  normalizing 
noise  sample  ana  tne  secona-oraer  or  fourth-oraer  estimate  can  performance  be 
improved.  This  result  could  impact  current  sonar  processing  methods. 
Simulations  were  presented  which  support  these  results. 
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APPENDIX  A 

SECOND-ORDER  ESTIMATES 


Let  the  power  spectrum  density  estimate  ot  equation  (14)  be  normalized  by 
a  noise-only  second-order  estimate  as  follows: 

ptFp)  =  xj/x 2  =  CXIQ)/Hi]2}/ {^/ [x(q)/HQj2 j>  (Al) 

q=l  q=l 

Asymptotically,  P t Fp j  wii]  also  converge  to  a  Gaussian  process,  unaer  H^, 
according  to  Cramer's  [lOj  convergence  proof. 


Employing  the  data  model  of  equation  (12)  we  obtain  the  expected  value  of 
equation  ( Al j  as 


ECP(Fp; j  .  1  +  L  SNR. 

The  components  of  the  variance  of  equation  (Al)  are  as  follows: 

VARLXj]  =  (1/n)  [2N2  +  4L  SN  ♦  (3L  -  L2)  S2  ] 

VAR(X2]  =  (1/n)  L2N2] 

Fj  =  1/N 

F2  =  -  U/N)(l  +  L  SNR) 

C0V[X1  X2]  =  (1/n)  2N2 

Therefore,  the  variance  of  equation  (Al )  is  given  by 


VAR[P(Fp)]  .  (1/n)  a8 


A-l 
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where,  ag  =  +  |_ 2)  SNR2  +  4L  SNR. 

As  SNR  approaches  zero  VARi_P ( Fp ^ j  also  approaches  zero.  This  result 
makes  sense,  because  when  the  SNR  equals  zero  P(Fp)  equals  one,  ana 
therefore  its  variance  woula  equal  zero. 

From  equation  (13)  we  obtain  the  result 
APD  =  1  -  <*>[-  L  SNR  n /  ag] 

As  SNR  approaches  zero  the  APD  approaches  .5.  However,  to  avoia 
complications  in  the  limit  we  will  assume  that  SNR  approaches  zero  but  aoes 
not  exactly  reach  it  under  Hq,  This  can  be  stated  more  precisely  as 
follows:  As  >  Hq,  SNR  *  6,  for  6  <<  1,  and  APD  »  .5  +  D(6),  where 
D(fi)  »  0  as  *  »  0.  This  result  also  establishes  the  limiting  false  alarm 
probability  of  .5  for  equation  (Al).  Therefore,  the  false  alarm  probability 
cannot  be  controlled  at  a  desired  low  level.  On  the  other  hana,  the  APD  woula 
have  higher  levels  comparea  with  the  FDK  of  case  2  as  a  function  of  SNR.  This 
may  not  be  of  interest,  however,  due  to  the  high  false  alarm  probability. 

If  we  now  let  the  PSD  estimate  of  equation  (14;  be  normalized  by  an 
uncorrelatea  noise-only  second-order  estimate,  say  from  an  adjacent  frequency 
location,  which  woula  give  C0V(X1 ,x2 )  =  0,  then  the  mean  aria  variance 
would  be  as  fol lows: 

E[P(Fp>J  =  1  +  l-  SNR 


yAR[P(Fp) j  =  (4/n)  q3> 
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where,  Q3  *  [1  +  2  L  SNR  +  ^3L  +  L2)  SNR2/4 j . 

These  results  give  an  APD  of  the  form 

APD  =  1  -a)  -  L  SNRv'TW2/v'?]/Q3/2}  •  (A2) 

The  false  alarm  probability  can  now  be  maintained  at  a  desired  level. 

But,  for  small  SNR  levels,  i.e.,  SNR^  <<  1,  the  APD  of  equation  (A2j  will  be 
lower  than  the  corresponding  APD  levels  of  equation  (15)  for  the  same  false 
alarm  probability.  This  is  true  because  of  the  factor ^/n ll  in  equation  (A2). 
If  L  equals  one,  then  the  APD  of  equation  (15)  will  always  be  higher 
independently  of  SNR.  For  high  SNR  there  may  be  values  of  L  where  the  APD  of 
equation  (A2)  is  higher  but  the  factor /n/2  must  be  overcome  before  this  is 
true. 
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APPENDIX  B 

VARIANCE  OF  THE  FDK  ESTIMATE  OF  CASE  1 


Here  we  present  the  major  points  for  the  derivation  of  the  variance  based 
on  Cramer's  method  [10J  for  the  following  equation 


W  =  VX21  ’ 

n  a  n  o  2 

wnere  X,  =  ( 1/n)  53  X(q)  »  and.  %  =  x(fl)  3 

q=l  q=l 


From  section  II  the  form  of  the  variance  is  known  and  it  can  be  expressed 


VARCKX(F  )]  -  VARIX^Fj  ♦  VAR(X21)F2  +  ZF^  C0V[X1X21] 


Tne  covariance  plays  an  important  role  in  minimizing  the  total  variance 


of  equation  (Bl). 


The  partial  derivatives  F^  and  F2  of  equation  (Bl)  are  as  follows: 


Fi  =  aK1(Fp)/aX1  x  =  J  =  1/ 1(  1  / n ) {! X ( q ) 4  -  (X(q)2)2]  +  (X(q)2)2j 

X21  3  hi 

where  tne  overbar,  wnen  used,  will  represent  the  expected  value, 


:2  -  3Ki(FP)/j*2i  *lmJi  "  _x(q) /j(1/n)  [x(q)  ‘  (X(q)  }  +  (X(q)  ' 


X21  =  *21 


B-l 
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The  terms  VAR(X^)  and  VAR(X2^)  can  oe  expressed  in  the  form 

VAR(X^)  =  (1/n)  [X(q)S  -  (X(q)4)2]  , 

VAR(X21)  =  (4/n)  -  (X(q)V  -  (X(q)2)4]  , 

if  the  terms  whicn  converge  to  zero  faster  than  (1/n),  i.e.,  0(l/n2),  are 
negl ected. 

The  covariance  is  defined  as 

C0VLX1X21]  -  E(X1X21)  -  E(XX)  E(X21)  . 

If  Xj  and  X21  are  independent,  i.e.,  E(X1X21)  =  E(Xj)  E(X21)  then 

cov(x1x21)  =.  0. 

From  the  definition  of  X^  and  X21  we  obtain 

C0V(X1X21/H1)  =  E{(l/n)£  X(q)4[(l/n)£  X(q)2]2} 

A  1  q-1  q=l 

-E[(l /n)£  X(q)4]E{[(l/n)£  X(q)2]2}  , 

q-1  q=l 

where  we  have  indicated  that  the  denominator  contains  signal  and  noise  by 
in  the  covariance  expression.  We  do  this  to  distinguish  between  case  1  and 
case  2  of  the  FOK  whicn  contains  only  noise  in  the  denominator.  By  neglecting 
terms  of  order  0(l/n^)  we  obtain 
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C0¥(X1X21/H1)  -  (2/n) |  X(q)6  X(q)?  -  X(q)4  (X(q)2)2|. 

Substituting  all  the  terms  into  equatin  (Bl)  we  obtain 

,  _  _  _  -4 

VARCK^Fp)]  =  (1/n)  J  [X(q)S  -  (X(q) V][X(qr] 

+  4[X(q)4  (X(q)2)2  -  ( X ( q ) 2) 4] ( X ( q ) 4 ) 2  (X(q7^)-8 
-  4[X(q)6  X(q) 2  -  X(q)4  (X(q)2)2]  X(q)4(X(q)2)"6 j 


In  order  to  express  equation  (B2)  in  the  desired  parameterized  form  the 
data  model  of  equation  (12)  is  suostituted  into  each  term.  But  instead  of 
doing  tnis  for  each  term,  only  the  general  term  of  equation  (B2)  will  be 
evaluated.  The  other  terms  will  follow  in  a  similar  way.  Therefore,  we  shall 
express  VAR(X^)  in  parameterized  form.  The  most  general  term  is 


VAR(X1)  =  X(q)S  -  (X(q)4}2  (B3) 

From  equation  (7)  the  first  term  of  equation  (B3)  is 

X(q)8  =  (1/M4)  0p8  |  E (x ( i x )  X(i2)  X(i3)  X(i4)  X(i5)  X(i&)  X(i;)  X(ig|J 


wnere  we  have  defined  tne  operator 


OP 


»  M-l  M-l  M-l  M-l  M-l  M-l  M-l  M-l 

8)  ZEEEZEED  |C0S(Fpi1)  cos(Fpi2) 

i^-0  i2»0  i3*0  i4»0  ig=0  ig=0  i^=0  ig=0 

cos(Fpi3)  cos(Fpi4)  cos(Fpig)  cos(Fpi6)  cos(Fpi7)  cos(Fpig) 


B-3 
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and  wnere  eacn  X(i)  is  also  a  function  of  q,  i.e.,  X(i,q)  =  N(i,q)  +  M(q)S(i,q), 
Dut  tnere  is  no  need  to  express  tnat  fact  in  the  equation. 

Since  N(i,q)  and  S(i,q)  are  Gaussian  OPgj  |  can  oe  expanded  as  follows: 

(1/M4)  0P8j  ECXd^  X ( i 2 )  Xd3)  X(i4)  X(i5)  X(1g)  X(i?)  X(ig)J 

*  105  N4  +  L  420  SN3+  630  L  S2  N2  +  L  420  N  S3  +  L  105  S4 

wnicn  is  not  a  function  of  the  transform  size  M.  By  defining  OP^j  j  for  the 
second  term  and  proceeding  in  a  similar  way  equation  (B3)  reduces  to 


X(q)8  -  (X(q)4)2  =  96  N4  +  384  L  S  H3  +  (612L  -  36L2)  S2  N2 
+  (420L  -  36L2)  S3  SN  +  (105L  -  9L2)  S4 

Tne  other  terms  are  obtained  similarly.  These  results  are  then 
substituted  into  equation  (B2)  to  obtain  equation  U7),  which  is  also  not  a 
function  of  the  transform  size  M. 
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Figure  2.  Ambient  Noise  witn  randomly  occurring  signal 


PROBABILITY  OF  OCCURRENCE  (L)  OF  THE  SIGNAL 

:,gUre  3.  Asymptotic  Prooaoility  of  Detection  vs  ProoaDil  ity  of  Occurrence 
(L)  of  the  signal .  FDK  (Case  1)  and  PSD. 
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Figure  4.  Real  Arctic  Data:  Top  -  Power  Spectrum  Density  Estimate; 

3ottom  -  Frequency  Domain  Kurtosis  Estimate 
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PROBABILITY  OF  OCCURRENCE  (L)  OF  THE  SIGNAL 

Figure  6.  Asymptotic  Probabil i ty  of  Detection  vs  Probabil  ity  of  Occurrence 
U;  of  tne  Signal.  F UK  (Case  2),  PSD,  ana  Fourth-order. 
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Figure  8.  Asymptotic  Probability  of  Detect 
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Figure  9.  Asymptotic  ProDaoility  of  Detection  vs  Si  goal -to-Noise  ratio. 
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Figure  10.  Simulation  of  FDK  of  Case  1.  PSD,  Real  Kurtosis,  Imaginary 
Kurtosis,  and  Magnitude  of  Kurtosis. 
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Figure  11.  Simulation  of  FDK  of  Case  2.  PSD,  Real  Kurtosis,  Imaginary 
Kurtosis,  and  Magnitude  of  Kurtosis.  L  =  .02. 
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Figure  13.  Simul at  ion  of  FDK  of  Case  2.  PSD,  Real  Kurtosis,  Imaginary 
Kurtosis,  and  Magnitude  of  Kurtosis.  L  =  .08. 


Kurtosis,  Imaginary  Kurtosis,  and  Magnitude  of  Kurtosis. 

(a)  SNR  -  .1,  (0)  SNR  -  .16,  (c)  SNR  -  .25  and  (d)  SNR  -  .4 
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